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Abstract 

We develop the theory of invariant random fields in vector bundles. The 
spectral decomposition of an invariant random field in a homogeneous vec- 
tor bundle generated by an induced representation of a compact connected 
Lie group G is obtained. We discuss an application to the theory of cosmic 
microwave background, where G = 50(3). A theorem about equivalence of 
two different groups of assumptions in cosmological theories is proved. 



1 Introduction 

This paper is inspired by Geller and Marinucci (2008). After reading the above 
paper and several physical books and papers cited below, the author realised that 
cosmological applications require the theory of random fields in vector bundles. 
A variant of such a theory is developed in Section 2, while an application to cos- 
mology is described in Section 3. 

According to vast majority of modern cosmological theories, our Universe 
started in a "Big Bang" . This term refers to the idea that the Universe has expanded 
from a hot and dense initial condition at some finite time in the past, and continues 
to expand now. 

*This work is supported by the Swedish Institute grant SI-01424/2007. 

^Division of Applied Mathematics, School of Education, Culture and Com- 
munication, Malardalen University, SE 721 23 Vasteras, Sweden. E-mail: 

anatoliy . malyarenko@mdh . se 



1 



As the Universe expanded, both the plasma and the radiation grew cooler. 
When the Universe cooled enough, it became transparent. The photons that were 
around at that time are observable now as the relic radiation. Their glow is 
strongest in the microwave region of the radio spectrum, hence another name cos- 
mic microwave background radiation, or just CMB. 

In cosmological models, it is usually assumed that the CMB is a single realisa- 
tion of a random field. A CMB detector measures an electric field E perpendicular 
to the direction of observation (or line of sight) n. Mathematically, n is a point 
on the sphere S 2 . The vector E(n) lies in the tangent plane, T n S 2 . In other words, 
E(n) is a section of the tangent bundle t, = (TS 2 , n,S 2 ) with 

7f(n,x)=n, nGS 2 , xeT n S 2 . 

It follows that cosmology uses the theory of random fields in vector bundles. 
A short introduction to vector bundles may be found in Geller et al (2009). It 
is not difficult to give a formal definition of a random field in a vector bundle. 
Indeed, let K be either the field of real numbers R or the field of complex numbers 
C. Let t, = (<f , n, T) be a finite-dimensional K-vector bundle over a Hausdorff 
topological space T . 

Definition 1. A vector random field on is a collection of random vectors { X(f ) : t G 
T} satisfying X(f) G n~ l {t), t G T. 

In other words, a vector random field on the base T of the vector bundle % is a 
random section of £ . 

To define a second order vector random field, assume that every space n~ l (t) 
carries an inner product. 

Definition 2. A vector random field X(t) is second order if E||X(r)||^_ 1( . f - ) < °o, 
tET. 

Next, we try do define a mean square continuous random field. The naive 
approach 

limE||X(s)-X(?)|| 2 = 

S—¥t 

does not work. If s, t G T with s ^ t, then X(s) and X(f) lie in different spaces. 
Therefore, the expression X(s) — X(f) is not defined. 

To overcome this difficulty, we extend an idea of Kolmogorov formulated by 
him for the case of a trivial vector bundle and published by Rozanov (1958) and 
Yaglom (1961). We start Subsection 2.1 by defining a scalar random field on the 
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total space S, which we call the field associated to the vector random field X(f). 
Then, we call X(Y) mean square continuous if the associated scalar random field 
is mean square continuous. 

Let G be a topological group acting continuously from the left on the base T . 
We would like to call a vector random field X(f ) wide sense left G-invariant, if the 
associated scalar random field is wide sense left G-invariant with respect to some 
left continuous action of G on the total space <§ . However, in general there exist no 
natural continuous left action of G on S '. In Definition 5, we define an action of G 
on S 1 associated to its action on the base space T . Then, we call a vector random 
field X(f ) wide sense left G-invariant, if the associated scalar random field is wide 
sense left G-invariant with respect to the associated action. 

In Subsection 2.2, we consider an important example of an associated action: 
the so called homogeneous, or equivariant vector bundles. They are important for 
us by several reasons. 

On the one hand, they have a natural associated action of some topological 
group G. Moreover, the above action identifies the vector space fibers over any 
two points of the base space. Therefore, all random vectors of a random field 
X(f ) in a homogeneous vector bundle lie in the same space. We prove that for 
homogeneous vector bundles, our definitions of mean square continuous field and 
invariant field are equivalent to usual definitions (3) and (4). 

On the other hand, the space of the square integrable sections of a homoge- 
neous vector bundle carries the so called induced representation of the group G. 
Therefore, we can use the well-developed theory of induced representations to 
obtain spectral decompositions of invariant random fields in homogeneous vector 
bundles. For an introduction to induced representations, see Barut and Raczka 
(1986). 

In Subsection 2.3 we consider mean square continuous random fields in ho- 
mogeneous vector bundles over a homogeneous space T = G/K of a compact 
connected Lie group G. In Theorem 1, we prove the spectral decomposition of a 
random field in a homogeneous vector bundle of the representation of the group G 
induced by an irreducible representation of its subgroup K. Here, we first meet the 
system of functions w Yy m (/) defined by (9), which form the orthonormal basis in 
the space of the square integrable sections of a homogeneous vector bundle under 
consideration. The spectral decomposition in Theorems 1-3 is given in terms of 
the above functions. 

In Theorem 2, we find the restrictions under which the spectral decomposition 
of Theorem 1 describes a wide sense G-invariant random field. Finally, Theo- 
rem 3 is a generalisation of Theorem 2 to the case when the representation of the 
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group G is induced by a direct sum of finitely many irreducible representations of 
the subgroup K. 

In Section 3, we apply theoretical considerations of Section 2 to cosmological 
models. Subsection 3.1 is a short introduction to the deterministic model of the 
CMB for mathematicians. In particular, we discuss different choices of local co- 
ordinates in the tangent bundle £ = (TS 2 ,n,S 2 ), and fix our choice. We explain 
both the mathematical and physical sense of the Stokes parameters I, Q, U, and V. 
The material of this Subsection is based on Cabella and Kamionkowski (2005), 
Challinor (2004), Challinor (2009), Challinor and Peiris (2009), Durrer (2008), 
and Lin and Wandelt (2006). 

The probabilistic model of the CMB is introduced in Subsection 3.2. We de- 
fine the set of vector bundles £, s = (<o s ,7t,S 2 ), s G Z, where the representation 
of the rotation group G = 50(3) induced by the representation W(g a ) = e lsa of 
the subgroup K = 50(2) is realised. In particular, the absolute temperature of 
the CMB, r(n), is a single realisation of a mean square continuous strict sense 
isotropic (i.e., 50(3) -invariant) random field in E,q, while the complex polarisa- 
tion, (<2±it/)(n), is a single realisation of a mean square continuous strict sense 
isotropic random field in £±2- Because any second order strict sense isotropic ran- 
dom field is automatically wide sense isotropic, Theorem 2 immediately gives the 
spectral decomposition of the above random fields. In the case of the absolute tem- 
perature, the functions (9) become familiar spherical harmonics, Y^ m , while in the 
case of the complex polarisation they become spin-weighted spherical harmonics, 
±?Xim- This fact explains our notation, wrYy m (f). The expansion coefficients are 
uncorrelated random variables with finite variance, which does not depend on the 
index m. In physical terms, the variance as a function of the parameter £ is the 
power spectrum. 

While studying physical literature, we have found that there exist various def- 
initions of both ordinary and spin- weighted spherical harmonics. The choice of a 
definition is called the phase convention. In terms of the representation theory, the 
phase convention is the choice of a basis in the space of the group representation. 
We made an attempt to describe different phase conventions in order to help the 
mathematicians to read physical literature. We also describe different notations 
for power spectra. 

Following Zaldarriaga and Seljak (1997), we construct the random fields E(n) 
and B(n). The advantage of this fields over the complex polarisation fields (Q ± 
it/) (n) is that the former fields are scalar (i.e., live in <fjo), real- valued, and isotropic. 
Moreover, only T(n) and E(n) may be correlated, while two remaining pairs are 
always uncorrelated. 
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Our new result is Theorem 4. It states that the standard assumptions of cosmo- 
logical theories (the random fields r(n), E(n), and B(n) are jointly isotropic) is 
equivalent to the assumption that ((Q — iU)(n),T(n), (Q + iU)(n)) is an isotropic 
random field in <^_2 © £o © £2- 

We conclude by two short remarks concerning Gaussian cosmological theo- 
ries and an alternative description of the CMB in terms of the so called tensor 
spherical harmonics. 

Note that we do not consider questions connected with statistical analysis of 
the observation data of the recent and forthcoming experiments. For an introduc- 
tion to this field of research, see Geller et al (2009) and the references herein. 

I am grateful to Professor Domenico Marinucci for useful discussions on cos- 
mology. 

2 Random fields in vector bundles 
2.1 Definitions 

Let (£2,#, P) be a probability space and let X(f) = X(f, fit)) be a vector random 
field in a finite-dimensional K- vector bundle B, = ($ , n, T). 

Definition 3. Let X(t,x) be the scalar random field on the total space S, defined 
as 

X(t,x) = (x,X(t)) 7[ -i {t) , teT, xen-\t). 
We call X(t, x) the scalar random field associated to the vector random field X(t). 

The field X(t,x) has the following property: its restriction onto 7T — 1 (?) is lin- 
ear, i.e., for any x, y G n^ 1 (t), and for any a, /3 G K, 

X(f,ax + /3y) = aX(t,x)+ fiX(t,y) P-a.s. (1) 

Definition 4. A vector random field X(f ) is mean square continuous if the asso- 
ciated scalar random field X(?,x) is mean square continuous, i.e., if the map 

P), (t,x)^X(t,x) 

is continuous. 



5 



Let H be a finite-dimensional K- vector space with an inner product (•,•). For 
any xE H, let x* be the unique element of the conjugate space H* satisfying 

x*(y) = (y,x), yeH. 

The mean value 

M(0 = E|X(0] 

of the mean square continuous random field X(f) is the continuous section of the 
vector bundle t, , while its covariance operator 

R(s,t) = E[X[s)®X(t)*} 

is the continuous section of the vector bundle £ Cg> £*. 

Because the scalar random field X(x) has property (1), it can be left invariant 
with respect to the associated action, only if the restriction of the associated action 
onto any fiber 1 (f ) is a linear invertible operator acting between the fibers. 
Moreover, the associated action must map the fiber n~ l (t) onto the fiber n~ l (gt). 

Definition 5. Let £ = (<o,7C,T) be a vector bundle, and let G x T — > T be a con- 
tinuous left action of a topological group G on the base space T. A continuous left 
action G x $ — ► $ of G on the total space $ is called associated with the action 
G x T — > T, if its restriction on any fiber n~ l (t) ls an invertible linear operator 
acting from 7T — 1 (?) to K~ 1 (gt). 

We are ready to formulate the main definitions of Subsection 2.1. 

Definition 6. Let £ = 7t, T) be a vector bundle, let G x T — > T be a continuous 
left action of a topological group G on the base space T, and let G x <f — > $ be 
an associated action of G on the total space $ . A vector random field X(f ) on ^ is 
called wide sense left G-invariant if the associated scalar random field X(t,x) is 
wide sense left invariant with respect to the associated action G x § — > $ , i.e., for 
all g E G, for all s,teT, and for all x E K~ l (s),y E n~ l (t) we have 

E[X(gs,gx)] = E\X{s,x)] t 
E[X(gs,gx)X(gt,gy)} = E[X(s,x)X(t,y)}. 

Definition 7. Under conditions of Definition 6, a vector random field X(f ) on t, 
is called strict sense left G-invariant if the associated scalar random field X(t,x) 
is strict sense left invariant with respect to the associated action G x $ — > <g , i.e., 
all finite-dimensional distributions of the random field X(t,x) are invariant under 
the associated action. 
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It is easy to see that any mean square continuous strict sense invariant random 
field is wide sense invariant. On the other hand, any Gaussian wide sense invariant 
random field is strict sense invariant. 

2.2 An example of associated action 

Let G be a topological group, and let K be its closed subgroup. Let T be the 
homogeneous space G/K of left cosets goK, go E G. An element g EG acts on T 
by left multiplication: 

g K ggoK. (2) 

Let W be a representation of K on a finite-dimensional complex Hilbert space 
H. Consider the following action of K on the Cartesian product G x H: 

k(g,x) = (gk,W(k- l )x). 

Denote the quotient space of orbits of the above action by $w. The projection 

7t:g w ^T : n(g,x)=gK 

determines the homogeneous, or equivariant vector bundle t, = (<%, n, T). 
Let t = goK E T. It is trivial to check that the action 

g(goK,x) = (ggoK,x) 

is associated to the action (2). 

Moreover, let X(f) be a random field in t s . All random vectors X(Y) lie in the 
same space H. By definition, the associated scalar random field X(t,x) = (x,X(Y)) 
is mean square continuous if and only if 

lim E\X(s,y)-X(t,x)\ 2 = 0. 

(*,y)-Kf,x) 

Let {ei,e2, • • • ,edi m //} be a basis in H. Put y = x = ej. Then we have 

\imE\Xj(s)-Xj(t)\ 2 =0, 

which is equivalent to 

timE\Xj{s) -Xj(t)\ 2 = 0. 
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It follows that 



dim// 



limE||X(5)-X(0|| 2 = limE £ \Xj(s) - Xj{t 



dim// 

£ limEI^W-Z/OI^O. 

7=1 



Conversely, let 



limE||X(>)-X(>)|| 2 =0. 



(3) 



Then, for any j = 1, 2, . . . , dim//, 



< lim sup E|X/ (01 

dim// 

< £ limsupE|X ; (s)-X 7 (f)| 2 = 0, 



./=i s ^ 

2 



thus, lim^ f E\Xj(s) -Xj(t) \ z = 0. It follows that 



lira E|X(,,y)-X(;,x)| 2 = lira E 

(J,y)->(f,x) (s,y)^(t,x) 



dim// 



dim// 

<2 £ lim Ely^W-x^WI 
£i (*j)-KM0 



0. 



We proved that in the particular case of a vector random field in a homoge- 
neous vector bundle our definition of mean square continuity is equivalent to the 
usual definition (3). In the same way one can easily prove that our definition of 
a wide sense G-invariant field is equivalent to the following equalities: for all s, 
t ET, and for all g E G we have 



(4) 



E[XW] = E[XW], 
E[X(^)®X*(g;)] = E[X(*)®X*(0]. 

The first equation is equivalent to the following equality 

E[X(s)] = E[X(f)], s,teT, 

because G acts transitively on T. Thus, the mean value of a wide sense G-invariant 
random field on T is constant. 
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2.3 The spectral decomposition of a vector random field over a 
compact homogeneous space 

Let G be a compact topological group, and let K be its closed subgroup. Let 
T be the homogeneous space G/K. Let W be a representation of K on a finite- 
dimensional complex Hilbert space H, and let % = T) be the correspond- 
ing homogeneous vector bundle. Let G (resp. .fif) be the set of all equivalence 
classes of irreducible unitary representations of G (resp. K). For simplicity, as- 
sume that K is massive in G (Vilenkin, 1968). This means that for all V E G and 
for all W E K the multiplicity of W in the restriction of V onto K is either or 1. 

First, consider the case when W is an irreducible unitary representation of K. 
Let dg be the Haar measure on G with J G dg = 1. Let L 2 (G,H) be the set of all 
measurable functions f : G — > H such that 

/ ll%)H 2 ds<°° 

and 

%fc)=W(fc- 1 )%), g EG, k E K. (5) 

To each f e L 2 (G,H), we associate the map s: T — > <%: s(g^T) = (g,f(g)). The 
above association is an isomorphism between L 2 {G,H) and the space L 2 (<%) of 
the square integrable sections of the homogeneous vector bundle t, . This space 
can be considered as a space of "twisted" functions on the base space T . If W is 
the trivial representation of K in H = C, then we return back to the standard space 
L 2 (G). The representation 

[^(g)s](t)=s(g-h) 

is the representation of G induced from the representation W of the subgroup K. 
We need the following facts about induced representations. 

1. Frobenius reciprocity: the multiplicity of V E G in is equal to the multi- 
plicity of W in V. 

2. The representation induced from the direct sum W\ © W2 © • • • © W# is the 
direct sum of representations induced from W\, W2, Wff. 

Let Gk(W) be the set of all V E G whose restrictions onto K contain W (nec- 
essarily once, because K is massive in G). For any V E Gr-(W), let iy be the em- 
bedding of H into the space Hy of the representation V. Let />i/ be the orthogonal 
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projection from Hy onto H. By the result of Camporesi (2005), any f G L 2 (G 1 H) 
can be represented by the series 

f{g) = d^W £ dimV f p v V(g- l h)i v f(h)dh. 
The above series converges in strong topology of the Hilbert space L 2 (G,H), i.e., 
l|f|l W) = dmW ^ dimV J G (j G Pvy(s- l h)ivmdhM) dg. 



VeG K (w) 

Ji,e 2 ,. 

basis in Hy with 



Fix a basis {ei , e 2 , . • • , e dimH } of the space H. Let {e^ ,e%\..., e^2 Hy } be a 



i V ej = e% l<j<dimW. (6) 

for some p > 0. Let //(g) = (f(g),ey) be the coordinates of f(g). Equation (6) 
means that W acts in the linear span of the dimW basis vectors of Hy that are 
enumerated without lacunas. Then we have 

i v f(A) = (0,...,0,/i(A),...,/ dimW (A),0,...,0). 

Let V m ,„(g) = (V(g)em\G^) be the matrix elements of the representation V. 
Then 

dim V dim W 

(y{g- l h)i v i(h)) p+J = £ ^+,-(g) £ v miP+n {h)f n (h) 

m= 1 n= 1 

and 

^ dimVdimW /. 

^^ = ^7 I dimV £ £ / f„(h)V m , p+n (h)dhV m , p+J (g). 

W VeG K (W) m=l n=l JG 

From now, let G be a connected compact Lie group, and let p : G — ► T denote 
a natural projection: p(g) = g£. Let be an open dense subset in G, and let 
(D G , J(g)) with 

J(g) = (0! (g) , . . . , dimG (g)) : Z> G - M dimG 

be a chart of the atlas of the manifold G with the following property: if k G K and 
both g and fcg lie in D G , then 0/(fcg) = 0j(g) for 1 < j < dimT. Then, (D T ,l{t)) 
with 

I(O = (0i(O,...,0d 1 mrW):^r->M dimr 
10 



is a chart of the atlas of the manifold T, and the domain Dj of this chart is dense 
in T. Let t £ Dj has local coordinates (6\,..., Odimr) m the chart (7). Then, the 
representation of the section s £ L 2 (<%) associated to f £ L 2 (G,H) has the form 

j dimVdimW n 

S i® = x^w E dimV E E Sn(t)V miP+n {t)dtV m , p+j (t) (8) 

VeG^(W) «=1 n=l yr 

where J? is the G-invariant measure on T with f T dt= 1, and 

Vm,p+«(0 = V m , p + n (0i, ... , OdimT, ^dimT+l' ' - ' ' ^dimG) - 

Introduce the following notation: 



dim V 

wYy m (f) = d di m (^",P+! ( f ) ' ^'",P+2 (0 > • • • i Vm,p+dhnw(t) ) ■ (9) 

Note that the correct notation must be w^Wm(t)> because functions (9) depend on 
the choice of a chart. In what follows, we use only chart (7) and suppress symbol 
I for notational simplicity. 

Equation (8) means that the functions { wYv m (t) '■ V £ Gx(W), 1 < m < dimV } 
form a basis in L 2 (<%), i.e., 

dimVdimW n 

S &) = E E E / s n(t)(wYvm)n{t)dt( W Y Vm )j{t) (10) 

Let X(t) be a mean square continuous random field in E, . Consider the follow- 
ing random variables: 

zl ] = J X n (t)( w Y Vm )n(t)dt, (11) 

where V £ Gx(W), 1 < m < dimV, and 1 < n < dim W. This integral has to be un- 
derstood as a Bochner integral of a function taking values in the space L 2 ^ (£2, 5, P) • 

Theorem 1. Le£ G be a connected compact Lie group, let K be its massive sub- 
group, let W be an irreducible unitary representation of the group K, and let t, be 
the corresponding homogeneous vector bundle. In the chart (7), a mean square 
continuous random field X(f) in t, has the form 

dim V dim W 

X j( t )= E E Y,z82(wYvm)M < 12 ) 

VeG K (W) m=l n=l 

where random variables Z^J have the form (11). 
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Proof. Let M(f) be the mean value of the random field X(f), and let R(ti,t 2 ) be 
its covariance operator. Denote the right hand side of (12) by Zj(t). We need to 
prove that 

EpE(0]=M(0 

and 

£[Z{h)®Z*{t 2 )]=R{t h t 2 ). 

Using (1 1), we obtain 



It follows that 

E[Zj(t)] = E 



E[z£?] = J E[X n (t)](wY Vm )n(t)dt. 



dimVdimW 

E E E zSiwYvmhit) 
VeG K (W) m=\ n=l 



dim V dim W 

III E\zS2]( w Y Vm )j[t) 

VeG K {W) '"=! » =1 

dim V dimW /. 

E EE/ E[Xn{t)}{wYvm)n(t)<lt{ W Yvm)j(t) 

E[*/(*)] 



(13) 



by (10). 

Similarly, 



P lz (v) 7 (v>) 



z iv] = JJ TxT R ( t ut2)(wYv,n)n(h)(wYv' m ')n'(h)dt 1 dt 2 . 



It follows that 



dim V dim V' dimly 



E[Z j (t l )Z f (t 2 )}= £ E E E E[Z^Z^J]( w 7v m ) J (^i)( l yW)/(?2) 

VjV'e^jffW) m= l »l'=l n,n'—l 

dim V dim V' dim W /. n 

= E E E E // m,t2)(wY Vm )n(ti) 

V,V'EG K {W) ™=\ m'=l n,n'=l JJTxT 

X (wYvi m i) n i{t 2 )At\At 2 (wYy m )j(t\)(wYyi m i)f{t 2 ) 

= Rjf(h,t 2 ). 

(14) 

□ 
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Denote by Vq the trivial irreducible representation of the group G. 

Theorem 2. Under conditions of Theorem 1, the following statements are equiv- 
alent. 

1. X(Y) is a mean square continuous wide sense invariant random field in 

2. X(f) has the form (12), where z£2, V G G K (W), 1 < m < dimV, \<n< 
dimW are random variables satisfying the following conditions. 

• IfV^Vo, thenE[Z { Z}}=0. 

• E[zS2z%J] = b\ v i8 mml R^J, with 

£ dimVtr[i?W]<oo. (15 ) 

Ve6 K (w) 

Proof. Let X(f) be a mean square continuous wide sense invariant random field 
in By Theorem 1, X(f) has the form (12). By (13), we have 

dim V dim W 

EP0«]= III E[zffi]( w F Vm ),.(0- 
Let g EG. Substitute gt in place of f to the last display. We obtain 

dim V dim W 

dim 1/ dim W dimV 

= E E E E t z -] E iU^faJWMO 

dim V dim W dim V 

= E E E E V em{g)E[zZ }(w Y Vm)j{t)- 
VeG K (W) m=i n=l 1=1 

The left hand sides of the two last displays are equal. Therefore, the coefficients 
of the expansions must be equal. 



£ Vt m (g)E[zjr>] = E[Z%>] 
i=i 

13 



Denote = (E[z£?] 5 . . . , EfZ^J). Then 

V + (g)M n V) =M n V \ gGG, 



where = V(g l ) T is the representation, dual to the representation V. It 

follows that either M. n = or the one-dimensional subspace generated by M.„ 
is an invariant subspace of the irreducible representation V + . In the latter case, 

V + must be one-dimensional. If V + is trivial, then V is also trivial, and M„ is 
any complex number. If V + is not trivial, so is V. Then, there exist g E G with 

V(g) + 1. It follows that Mf } = V(g)M n V) = 0. 
By (14), we have 

dim V dim V' dimW 



Rjf(h,t 2 )= £ E E E E[z^z^J]( w y Vm ),(?i)( l vW)/fe). 

V,VeGA-(lV) m= l m!=\ n,n' = l 

It follows that 

dim V dim V' dim IV 



R jf (gti,gt 2 )= £ E E E EfZ^Z^jj^F^),-^!)^^)/^) 

V.V'eG^lV) m= l m '=l n,n'=l 

dimVdimV' dimVK 

= E LEE ^7^nn^ m i n i\ 
V,V'eG/s:(lV) m= l m'=l n,n'=l 

dimV dimV" 

x £ E v L'i'{s){wYvi')f{h) 

t=\ £>=1 

dimVdimV' dimW dim V dim V' 



= E E E E E E VntisMM 

V,V'eGi((W) m=l m' = l n,n'=l £=l £'=l 

x EiZ^Z^iwYv^jih^wYv^jitj). 
By equating the coefficients of the two expansions, we obtain 

dim V dim V' 



E E ^)^(y)E[^2^] = E[zS?z23] 

l=\ £' = \ 



Let ^ be the matrix with elements 

nn 
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Then, 

(vW)( g )pjrWir> S ea 

It follows that either pf^J V ' is zero matrix or the one-dimensional subspace gen- 
erated by P^f ^ is an invariant subspace of the representation V + <8> V'. In the 
latter case, the representation V + <8> V' contains an one-dimensional irreducible 
component, say Y. If V is trivial, then V = V' and "V acts in the one-dimensional 
subspace generated by the unit matrix {P\f ^)mmf- If 1S not trivial, there ex- 
ist g G G with Y(g) ^ 1. It follows that P [ J/' ] = r[g)P^f\ so P^P is zero 
matrix. So, 



E[Z^ Z^J] — 8y V i8 mm iR^J . 

Let ?o G T be the left coset of the unit element of G. We may assume to G Dj 
(otherwise use a chart {gDT : \{g^ l t)) for a suitable g G G). Then 



and 



(u-Yvn, )./('(>) = \/ 



dim V dim IV 

*/(«>) = ELI 32?(wiv»MO 



1 



dim IV 



£ Vdimy£z^. 



VdimVF 



veG K (w) 



n=\ 



It follows that 



and 



E| ^ (?0)|2 = dim^ E dimydim W E|Z^| 2 

= £ dimV*M 



dim W 



£ dimytr[i?( v )]= £ E|X / -(?o)| 2 <°°. 

VeG K (W) J= l 

Conversely, let Z^n , V G Gk(W), 1 < m < dimV, 1 < n < dimW be random 
variables satisfying conditions of Theorem 2. Consider random field (12). Then, 
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its mean value is 

r(V 

E[Xj(t 



E[zjn, V eG K (W), 
0, otherwise, 



which is constant. Note that Vq G Gk(W) if and only if W is trivial (by Frobenius 
reciprocity). 

The correlation operator of the random field (12) is 



dim V dim V' dimW 



Rjf{h,h)= 52 H H H E[ZmnZ m/n ,](wYvm)j(ti)(wYvi m ')f(t2) 
V,V'<eG k {W) m=l m'=\ n,n'=l 



dimW , , dimV 



= ]£ H R nn' H {w Y Vm)j{h){w Y Vm)j'{h) 
VeG K (W)n,n'=l m=l 

j dimly 

= dhnw I dimV I Ri >P + j,P + As^g2), 

VeG K (W) n,n'=\ 

where gi (resp. g2) is an arbitrary element from the left coset corresponding to 
t{ (resp. ti). The terms of this functional series are bounded by the terms of 
the convergent series (15), because \V p+ j^ p+ f(g\ l g2)\ < 1. Therefore, the series 
converges uniformly, and its sum is continuous function. This means that X(t) is 
mean square continuous. 
For any g e G, we have 

i dimly 
Qlm W VeG K (W) n,n'=l 

= Rjf(h,t2), 

so X(t) is invariant. □ 

Assume that W is not necessarily irreducible representation of K in a finite- 
dimensional complex Hilbert space H. Because K is compact, the representa- 
tion W is equivalent to a direct sum W\ © W2 © • • • © Wn of irreducible unitary 
representations of K. The representation induced by W is a direct sum of rep- 
resentations induced by W^, 1 < k < N. It is realised in a homogeneous vector 
bundle £ = £1 © £2 © • • • © where t, k is the homogeneous vector bundle that 
carries the irreducible component Wk- 
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Let X(t) be an invariant random field in % . Denote the components of X(f ) by 

Xj (t), 1 < k < N, 1 < j < dimWfc. Denote by Pk the orthogonal projection from 
H onto the space Hk where the irreducible component Wk acts. 

Theorem 3. Under conditions of Theorem 1, the following statements are equiv- 
alent. 

1. X(f ) is a mean square continuous wide sense invariant random field in t,. 

2. X(f ) has the form 

dimKdimWj. 

x j \')= E E E Z^ k \w k Yv m )j(t), (16) 

VeG K {W k ) m = l n=\ 

where Z { J n k \ 1 < k < N, V e G K (W k ), 1 < m < dimV, 1 < n < dimW k are 
random variables satisfying the following conditions. 

• IfV^Vo, thenE[2% n k) }=0. 

• E[zi?¥^f } ] = ^v'5 m «'4l,k" with 

N 

£ £ dimVtr[jy?( v )j\]<oo. 

Proof. Use mathematical induction. The induction base, when N = 1, is Theo- 
rem 2. Assume the induction hypotheses: Theorem 3 is proved up to N — 1 . 

Let X(t) be a mean square continuous wide sense invariant random field in t, . 
Then the field 

Yl (t) = {x® (t),... ,^ (0, . . (0, • • • (0), 

is a mean square continuous wide sense invariant random field in £i © • • • © <^v-i • 
By the induction hypotheses, 

dimVdimW^ 

Xj\t)= E E E z ™ kiym)X0> i<*<jv-i, 
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where E[z£?] = unless V ^V and E[Z$p2%j '] = ^vv' 8 mm' R krflri^ witn 



AV,N-l). 



£ £ dimVtr^J?^- 1 ^] < 

The field 



y 2 (0 = (< ) (0,.-.,^W) 

is a mean square continuous wide sense invariant random field in §N. By Theo- 
rem 2, 

dim 1/ dim 1% 

xf\t) = III zt N \w N Yv m )j{t)i 



V<EG K (W N ) 



m=\ n=l 



where E[Z%P] = unless V ^V and EfZjS^Z^J^] = 8 VV > 8 mrr jR^J' , with 

£ dimVtr[tf<™>] < 

veG K (w N ) 

The matrix R^nk'ri w * m e l ements 



kn.k'n L In 1m' J 

obviously satisfies conditions of the second item of Theorem 3. 

Conversely, let Z^ n k) , \<k<N,Ve G K (W k ), 1 < m < dimV, \<n< 6imW k 
be random variables satisfying conditions of Theorem 3. Consider random field 
(16). Its mean value is obviously constant. Its correlation operator is 

dimWfcdim% dimV 

Rff\tut 2 )= £ £ £ R ( k V n \, n , £ (w k Yv m )j(h)(w kf Yv m )f(t 2 ) 

VeG K {W k )nG K {W k ,) n=\ n'=\ m=\ 



y dimWj-dimWj,/ 

v/dimlV t dimHV . „"-.,.,,, , ~ ~, 



V€Gr(Wt)n6 t (HV) 



with the same notation as in proof of Theorem 2. The uniform convergence of the 
above series and the invariance of the field (16) is proved exactly in the same way 
as in proof of Theorem 2. □ 
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3 Application to cosmology 



3.1 The cosmic microwave background 

Let E(n) e T n S 2 be the electric field of the cosmic microwave background. From 
the observations, we define the intensity tensor. In physical terms, the intensity 
tensor is 

^ = C(E(n)<g)E*(n)), 

where (•) denote time average over the historical accidents that produced a par- 
ticular pattern of fluctuations. Assuming ergodicity, time average is equal to the 
space average, i.e., average over the possible positions from which the radiation 
could be observed. The constant C is chosen so that & is measured in brightness 
temperature units (in these units, the intensity tensor is independent of radiation 
frequency). It will be ignored in what follows. 

Introduce a basis in each tangent plane T n S 2 . Realise S 2 as { (x,y, z)el 3 :i 2 + 
y 2 + z 2 = 1} and define the chart (C//,h/) as Ui = S 2 \ {(0,0, 1), (0,0, -1)} and 
h/(n) = (0(n), <p(n)) e R 2 , the spherical coordinates. Let 50(3) be the rotation 
group in IR 3 . For any rotation g, define the chart (U g ,h g ) as 

U g = gUi, h g (n)=h I (g- 1 n). 

The sphere S 2 , equipped with the atlas { (U g ,h g ) : g £ 50(3) }, becomes the real- 
analytic manifold. The local 0-axis in each tangent plane is along the direction of 
decreasing the inclination 0: 

d 

The local <p-axis is along the direction of increasing the azimuth q>: 

e«p = (l/sin0)^. 

With this convention, eg, e<p, and the direction of radiation propagation — n form a 
right-handed basis. This convention is in accordance with the International Astro- 
nomic Union standard. The orthonormal basis (ee,e<p) turns S 2 into a Riemannian 
manifold and each tangent plane T n S 2 can be identified with the space M 2 . 

In the just introduced basis, the intensity tensor becomes the intensity matrix: 

& ab = (E a (n) ®E£(n)), a,b e {0,<p}. 
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The rotations about the line of sight together with parity transformation n — > — n 
generate the group 0(2) of orthogonal matrices in IR 2 . The action of 0(2) on 
the intensity matrix extends to the representation g i— > g£/g~ l of 0(2) in the real 
4-dimensional space of Hermitian 2x2 matrices srf with inner product 

This representation is reducible and may be decomposed into the direct sum of 
three irreducible representations. 

The standard choice of an orthonormal basis in the spaces of the irreducible 
components is as follows. The space of the first irreducible component is gener- 
ated by the matrix 

1 1/1 0\ 
2 a ° = 2> l) 

The representation in this space is the trivial representation of the group 0(2). 
Physicists call the elements of this space scalars. 

The space of the second irreducible component is generated by the matrices 

1 1/0 1\ 1 1/1 \ 

2 ai = 2\l 0j' 2 a3 = 2v0 -I)' 

Let g a e 50(2) with 



It is easy to check that 

ga(?ig~a l = cos(2a)t7i + sin(2a)cr 3 , 
ga03g a = -sin(2a)ai +cos(2a)c7 3 . 

The elements of this space are symmetric trace-free tensors. 

Finally, the space of the third irreducible component is generated by the matrix 

1 1 /0 -i\ 
2° 2= 2\i oj 

The representation in this space is the representation g i— > detg of the group 0(2). 
Physicists call the elements of this space pseudo-scalars (they do not change under 
rotation but change sign under reflection). The matrices <j\, 02, and (73 are known 
as Pauli matrices. 



ga 



cos a sma 
-sin a cos a 



(17) 



20 



The standard physical notation for the components of the intensity matrix in 
the above basis is as follows: 



& = l - (7(J + Uoi + Va 2 + Qo 3 ), 

or 

1 / I + Q U-iV 



The real numbers /, Q, U, and V are called Stokes parameters. Their physical 
sense is as follows. / is the total intensity of the radiation (which is directly propor- 
tional to the fourth power of the absolute temperature T by the Stefan-Boltzmann 
law). On the tangent plane T n S 2 , the tip of the electric vector E(n) traces out an 
ellipse as a function of time. The parameters U and Q measure the orientation of 
the above ellipse relative to the local 0-axis, eg. The polarisation angle between 
the major axis of the ellipse and eg is 

1 _j U 

* = 2 tan 0' 

and the length of the major semi-axis is (Q 2 + t/ 2 ) 1 / 2 . The last parameter, V, 
measures circular polarisation. 

According to modern cosmological theories, the polarisation of the CMB was 
introduced while scattering off the photons by charged particles. This process can- 
not induce circular polarisation in the scattered light. Therefore, in what follows 
we put V = 0. 

The physics of the CMB polarisation is described in Cabella and Kamionkowski 
(2005), Challinor (2004), Challinor (2009), Challinor and Peiris (2009), Durrer 
(2008), Lin and Wandelt (2006), among others. Of these, Challinor and Peiris use 
the right-hand basis, while the remaining authors use the left-hand basis, in which 
ee = d/d0. In what follows, we use the left-hand basis e#, e<p, — n with 

ee = ^, e, = (l/sin0)A. (19) 



3.2 The probabilistic model of the CMB 

The absolute temperature, T(n), is a section of the homogeneous vector bundle 
<^o = (<£b, 7t, S 2 ), where the representation of the rotation group G — 50(3) induced 
by the representation W(g a ) = 1 of the massive subgroup K = 50(2) is realised. 
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The representations V of the group G are enumerated by nonnegative integers 

£ = 0, 1 , The restriction of the representation Vp onto K is the direct sum of the 

representations e lma , m = —£, —£+!,...,£. Therefore we have dimVf = 21 + 1 
and \m\ < I. By Frobenius reciprocity, Gk(W) = {Vo, Vi, . . . ,VJ?, . . . }. 

The representations e lma of ^ act in one-dimensional complex spaces H m . To 
define a basis in the space of the representation Vg, choose a unit vector e m in 
each space H m . Each vector e m of a basis can be multiplied by a phase e 1C(m . The 
choice of a phase is called the phase convention. 

Any rotation g e 50(3) is defined by the Euler angles g = (<p, 0, y/) with <p, 
i// G [0,2tt] and e [0,^]. The order in which the angles are given and the axes 
about which they are applied are not subject of a standard. We adopt the so called 
zxz convention: the first rotation is about the z-axis by i//, the second rotation 
is about the jc-axis by 6, and the third rotation is about z-axis by (p. Note that 
the chart defined by the Euler angles satisfies our condition: the first two local 
coordinates (<p, 0) are spherical coordinates in S 2 (up to order) with dense domain 



The matrix elements of the representation Ve are traditionally denoted by 



and called Wigner D-functions. The explicit formula for the Wigner D-f unction 
depends on the phase convention. Choose the basis {e ra : — £ < m < £} in every 
space to obtain 



U,. 



D y J,((p,6,iif) = (Ve((p,0,\i/)e m ,e n ) H (e) 




J e )(a)-( iyn / (l + m)!(l-m]j ■ 2 e ( a /?) 




rn , i V r J \r — m — n 




cot 



2r—m—n 



(0/2). 



(20) 



The following symmetry relation follows. 




(0) = (-l)"-"^(0) 



(21) 
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In this particular case, formula (9) takes the form 



The functions in the left hand side form an orthonormal basis in the space of the 
square integrable functions on 5 2 with respect to the probabilistic 50(3) -invariant 
measure. It is conventional to form a basis with respect to the Lebesgue measure 
induced by the embedding 5 2 C M? which is An times the probabilistic invariant 
measure, and omit the first subscript: 

This is formula (A4.40) from Durrer (2008) defining the spherical harmonics. 

In cosmological models, one assumes that T(n) is a single realisation of the 
mean square continuous strict sense 50(3) -invariant random field in the homo- 
geneous vector bundle £ = {<S'q,tz,S 2 ). It is custom to use the term "isotropic" 
instead of "50(3) -invariant". By Theorem 2, we have 

l 

T ( n ) = E E z te y £m (n), 

£=0m=-£ 

where E[Z £m ] = unless I = and E[Z^„Z^ m /] = dui8 mm iR^ with 

£(2£+l)J?^ <oo. 

This formula goes back to Obukhov (1947). 

Physicists call Zg m s the expansion coefficients, and Ry) the power spectrum of 
the CMB. Different notations for the expansion coefficients and power spectrum 
may be found in the literature. Some of them are shown in Table 1 . 

In what follows, we use notation of Lin ans Wandelt (2006). In this notation, 
the expansion for the temperature has the form 

e 

T(n) = £ £ a T/m Y em (n)- (22) 

Since r(n) is real, the coefficients aji, n must satisfy the reality condition which 
depends on the phase convention. For our current convention, when the Wigner 
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7c 

*^£m 




Cabella and Kamionkowski (2005) 


a £m 




Challinor (2005), 






Challinor and Peiris (2009) 


Tim 




Durrer (2008), 






Weinberg (2008) 


&£m 


Q 


Lin and Wandelt (2006), 






Zaldarriaga and Seljak (1997) 


ar,£m 


Ct£ 


Kamionkowski et al (1997) 


a £m 


r T 



Table 1: Examples of different notation for temperature expansion coefficients 
and power spectrum. 

J-function is determined by (20), we have 

= ]j^e-^ (-1)^1(6) 

= {-i) m Y^{e^). 

Here we used the symmetry relation (21). The reality condition is 

ar,l -m = ( — 1 ) m <^T,im- (23) 

This form of the reality condition is used by Cabella and Kamionkowski (2005), 
Challinor (2004, 2009), Challinor and Peiris (2009), Durrer (2008), Kamionkowski 
et al (1997), Lin and Wandelt (2006), among others. 

Introduce the following notation: m + = max{m,0} = (|m| +m)/2, mT = 
max{— m,0} = (|m| — m)/2. We have (— m)~ = m + and m + — m = mT . If we 
choose another basis, { (— l) m e ra : — i < m < £}, then the Wigner D-function, 
Dml((p,0, y/), is multiplying by ( — 1)'" + " , and we obtain 

Y^ m (e,cp) = /^(-i) ( - mr e-"^, (e) 

= \/^(- 1 )" 7+ ^ e " m,P (- 1 ) md io(e) 
= Y em (6,<p). 
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The modified reality condition is 



o~T,l-m — a T,lm- (24) 

This form of reality condition is used by Geller and Marinucci (2008), Weinberg 
(2008), Zaldarriaga and Seljak (1997), among others. 

Let Tq = E[r(n)]. The temperature fluctuation, Ar(n) = T(n) — Tq, expands 

as 

c° £ 

AT(n) = Y, H a T ,£mYe m (n). 

e=im=-e 

The part of this sum corresponding to £ = 1 is called a dipole. When analysing 
data, the dipole is usually removed since it linearly depends on the velocity of the 
observer's motion relative to the surface of last scattering. 

The complex polarisation is defined as Q + iU. It follows easily from (18) 
that any rotation (17) maps Q + \U to e 2ia (Q + iU). Then, by (5), (g + iZ7)(n) 
is a section of the homogeneous vector bundle <^_2 = (^-2,^,5 2 ), where the 
representation of the rotation group G = 50(3) induced by the representation 
W(g a ) = e~ 2ia of the massive subgroup K = 50(2) is realised. By Frobenius 
reciprocity, G K (W) = {V 2 ,V 3 , ...,V £ ,...}. 

In general, let jgZ, and let = {<S'- s ,Tl,S 2 ) be the homogeneous vector 
bundle where the representation of the rotation group 50(3) induced by the rep- 
resentation W(g a ) = e~ lsa of the massive subgroup 50(2) is realised. In the 
physical literature, the sections of these bundle are called 

• quantities of spin s by Challinor (2009), Challinor and Peiris (2009), Geller 
and Marinucci (2008), Newman and Penrose (1966), Weinberg (2008) among 
others; 

• quantities of spin —s by Cabella and Kamionkowski (2005), Lin and Wen- 
delt (2006), Zaldarriaga and Seljak (1997) among others; 

• quantities of spin \s\ and helicity s by Durrer (2008) among others. 

Let g = (0, <p, y) be the Euler angles in 50(3). Put Y = 0. Then, n = (0, <p,0) 
are spherical coordinates in 5 2 . By (10) we obtain 

(Q + i£/)(n) = £ £ a- 2l tm-2Yem[n), 

£=2m=-£ 
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where 

a-2,im= (2 + it/)(n)- 2 ^m(n)dn 

Js 2 

and, by (9), 

-2Y im (e,(p) = V2£+lD^_ 2 ((p,e,0). 

The functions in the left hand side form an orthonormal basis in the space of the 
square integrable sections of the homogeneous vector bundle <^_2 with respect to 
the probabilistic 50(3)-invariant measure. 

There exist different conventions. The first convention is used by Durrer 
(2008) among others. In this convention, a basis is formed with respect to the 
Lebesgue measure induced by the embedding S 2 C M 3 which is 4n times the 
probabilistic invariant measure and the sign of the second index of the Wigner 
D-function is changed (because we would like to expand Q + iU with respect to 
lYlm)'- 

-2YU0,cp) = y/^^D^(<p,e,0). 

In the general case, for any s EZ, this convention reads (Durrer (2008), for- 
mula (A4.51)) 

s Y em (0,(p) = y^pI>W_,(9,e,0). (25) 

These functions are called spherical harmonics of spin s or the spin-weighted 
spherical harmonics. They appeared in Gelfand and Shapiro (1952) under the 
name generalised spherical harmonics. The current name goes back to Newman 
and Penrose (1966). Note that the spin- weighted spherical harmonics are defined 
for £ > \s\ and |m| < t. 

The second harmonic convention is used by Lin and Wandelt (2006), Newman 
and Penrose (1966), among others. It reads as 

s Y tm (e,cp) = (_i)»y^lz)g_(9,e,o). 

Both conventions are coherent with the following phase convention: 

sYlm = {-\) m+> ''- s Yi - m . (26) 

In particular, for s = we return back to the convention Y^ m = (— l) m 7| _ m corre- 
sponding to reality condition (23). 
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To produce the harmonic convention coherent with the phase convention 



corresponding to reality condition (24), one must multiply the right hand side of 
the convention equation by (— l) m . Thus, the modified first convention, used by 
Weinberg (2008) among others, is 



while the modified second convention, used by Geller and Marinucci (2008), 
among others, is 



s Yae,c P ) = (-i) m y 2 ^D^_ s (( P ,e,o). 

In what follows, we use the convention (25). The explicit expression for the 
spherical harmonics of spin s in the chart determined by spherical coordinates 
follows from (20) and (25): 



Y (Q<n\-( If / (2£+l)(£ + m)\(£-m)\ u ^ 



min{£+m,£— s} / g 



X 

r=max{0,m— s} 



(27) 



The decomposition of the complex polarisation takes the form 

oo £ 

(<2 + iC/)(n) = £ £ «Vm2^m(n), 

£=2m=-i 



where 



ai i i m = l(Q + iU)(n) 2 Y em (n)dn, 

J s 



is 2 

while the decomposition of the conjugate complex polarisation is 
(G-i£/)(n) = £ a -2,fm-2Ye m (n), 

e=2m=-£ 
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Source 


Expansion coefficients 


Durrer (2008) 


(±2) 
a lm 


Lin and Wandelt (2006), 
Zaldarriaga and Seljak (1997) 


a±2,£m 


Weinberg (2008) 


ap.lm 



Table 2: Examples of different notation for complex polarisation expansion coef- 
ficients. 

where 

a-2im= (<2-i^)(n)-2^m(n)dn. 

Js 2 

In cosmological models, one assumes that (Q + iU)(n) is a single realisation 
of the mean square continuous strict sense isotropic random field in the homoge- 
neous vector bundle % 2 - Isotropic random fields in vector bundles £ s , s E Z were 
defined by Geller and Marinucci (2008). By Theorem 2, we have 

00 £ 

(Q + i£/)(n) = £ £ a 2 ,tm2Yim(n), (28) 

e=2m=-£ 

where E[a 2 ,£ m ] = and E[a 2 /m«2/m'] = S U 'd mm 'C 2 £ with 

00 

£(2£+l)C 2f <°c 

1=2 

Different notations for the complex polarisation expansion coefficients a± 2 ,£m 
may be found in the literature. Some of them are shown in Table 2. 

In what follows we use the notation by Lin and Wandelt (2006). The expansion 
for the conjugate complex polarisation has the form 

00 e 

(Q-i£/)(n) = £ £ a-vm -2^«(n). (29) 

£=2m=-£ 

Since Q(n) and t/ (n) are real, the coefficients a 2 ^ m and a- 2 / m must satisfy 
the reality condition which depends on the phase convention. We agreed to use 
the first harmonic convention (25). Therefore, our phase convention is (26), and 
the reality condition is 

a~2,lm — (— l) m «2,/ -m- (30) 
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Along with the standard basis (19), it is useful to use the so called helicity 
basis. Again, there exist different names and conventions. Durrer (2008) defines 
the helicity basis as 

e+ = -^(ee-ie 9 ), e = -^=(e e +ie 9 ), 
while Weinberg (2008) uses the opposite definition 

e + = -^=(e fl +ie 9 ), e_ = -^(e e -ie 9 ). 
Challinor (2005) and Thorne (1980) use notation 

m= -^(ee + ie,,,), m* = -^=(e e -ie 9 ), 
while Challinor and Peiris (2009) use notation 

l / . N. l / . V 

m + = -=(e e +ie <?) ), m_ = --=(e -ie 9 ) 

and call these the null basis. We will use the definition and notation by Durrer 
(2008). 

The helicity basis is useful by the following reason. Let 3 be a covariant 
derivative in direction — v2e_: 



9 = V 



•V5e_- 



Let C°°(% s ) be the space of infinitely differentiate sections of the vector bundle 
£ s . Durrer (2008) proves that for any s f E C°°(£ s ) we have 

"^(""'-^-seiy- 

In particular, put s f = s Y^ m . Using (27), we obtain 

9 s Yi m = ^(£-s)(£ + s+l) s+ iY em . 

For s > and £ = j, the spherical harmonic s +iYe m is not defined and we use con- 
vention y/(£ -£)(!£ + \)i+\Y im = 0. Then, 9: C°°(% s ) -»■ C"(4+i). Therefore, 
9 is called the spin raising operator. Moreover, the last display shows that the 
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restriction of 9 onto the space £ > s, is an intertwining operator between 
equivalent representations Vi. 

The adjoint operator, 9*, is a covariant derivative in direction — v2e+: 



For any s / G C°°(^ s ) we have 



In particular, 



Q%Y £m = -^(£ + S )(£- S+ l) s _ 1 Y em . 



For s < and £ = —5, the spherical harmonic s -\^im is not defined and we use 
convention ^/(£-£)(2£ + = 0. Then, 9* : C°°(&) -> C°°(&_i). There- 

fore, 9* is called the spin lowering operator. Moreover, the last display shows 
that the restriction of S* onto the space £ > —s, is an intertwining operator 
between equivalent representations Vg. 

Zaldarriaga and Seljak (1997) introduced the following idea. Assume for a 
moment that 

f (2l+l)(l + 2)! ^ 

£ (/-2)! C ^ <0 °- (31) 

Then, it is possible to act twice with 9 on both hand sides of (29) and to inter- 
change differentiation and summation: 

9 2 (£-i£/)(n)=g 2 £ £ a_ 2 ,te-2^m(n) 

f=2m=-f 

= ^ ^ a_2/ m S 2 -2^m(n) 

f=2ra=-< 

f < / (1 + 2)! 

= 2- 2, y (£-2)! a_2,£m 

By the same argumentation, we have 



(r) 2 (£ + it/)(n) = f; £ J^^a 2/m Y im (n). 



e=2m=-e 
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Unlike complex polarisation, the new random fields are rotationally invariant and 
no ambiguities connected with rotations (18) arise. However, they have complex 
behaviour under parity transformation, because Q(n) and U (n) behave differently 
(Lin and Wandelt (2006)): Q has even parity: Q{— n) = Q(n) while U has odd 
parity: U(-n) = -U(n). 

Therefore, it is custom to group together quantities of the same parity: 

£(n) = -i((S*) 2 (e + it/)(n) + g 2 (2-i[/)(n)), 

5(n) = -i((5*) 2 (2 + it/)(n)-g 2 (e-i[/)(n)). 

The random fields E(n) and B(n) are scalar (spin 0), real- valued, and isotropic. To 
find their behaviour under parity transformation, follow Lin and Wandelt (2006). 
Notice that if n has spherical coordinates (6,<p), then — n has spherical coordi- 
nates 9' = % — and <p' = <p + %. Therefore, 

~do~'~ de' d<p>~d(p' 

Because (Q + iU)(-n) = (Q-iL/)(n), we obtain 

flWfl+ilOf-) = ( 2 co. 9 '- A + jL) (e + i[ /)(-„) 

= -3(e-it/)(n) 

and 

( ( 57) 2 (e + ^)(^ 

= d 2 (Q-iU)(n). 

Similarly, we have (&) 2 (Q + iU)(-n) = (9*) 2 (g - it/)(n). Therefore, 

£(-n) = -i((r) 2 (e+i[/)(-n)+g 2 (e-it/)(-n)) 
= -^((r) 2 (e-it/)(n)+g 2 (e+it/)(n)) 

= E(n) 
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and 

5(-n) = -i((g*) 2 (e + it/)(-n)-S 2 (e-it/)(-n)) 

= -^((r) 2 (e-it/)(n)-g 2 (e+it/)(n)) 

= -B(n). 

It means that E(n) has even parity like electric field, while B(n) has odd parity 
like magnetic field. 

The spectral representation of the fields E(n) and B(n) has the form 

oo e 

£=2m=-i 

oo I 

B(n) = £ J a Sem Y £m {n), 

e=2m=-£ 

where 



1 (£ + 2)1. 

a E,£m = ~2 y 7^jl fl 2,ta + fl-2,l 



(32) 



1 (£ + 2)\ 

a B,£m = ~2i\ j£~2yS a2/ ' m ~ a ^ 2/m ^' 

It is convenient to introduce the fields E(n) and B(n) as 

£=2m=-£ 

oo £ 

^=2ra=-f 

with 

a^ m = -^(a 2/m -a_ 2/m ). 

The random fields £(n) and 5(n) are scalar (spin 0), real- valued, and isotropic. 
Moreover, E(n) has even parity, while B(n) has odd parity. The advantage of 
£(11) and B(n) is that their definition does not use assumption (31). The expansion 



(33) 
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Source 


Fields 


Multipoles 


Challinor (2005), 






Challinor and Peiris (2009) 






Durrer (2008), 


£(n),@(n) 




Geller and Marinucci (2008) 


Ie, Im 




Lin and Wandelt (2006), 


E(n),B(n) 


a E,hm a Bfin 


Weinberg (2008), 






Zaldarriaga and Seljak (1997) 




a E,bm a Bfin 



Table 3: Examples of different notation for the fields E(n) and B(n) and its ex- 
pansion coefficients. 



coefficients ciE/ m are called electric multipoles, while the expansion coefficients 
a Bim are called magnetic multipoles. 

Different notations for the fields E(n) and B(n) and electric and magnetic 
multipoles may be found in the literature. Some of them are shown in Table 3. In 
what follows, we use notation by Lin and Wandelt (2006). 

We prove the following theorem. 

Theorem 4. Let T(n) be a real-valued random field defined by (22). Let (<2± 
it/)(n) be random fields defined by (28) and (29). Let E(n) and B(n) be random 
fields (32) whose expansion coefficients are determined by (33). The following 
statements are equivalent. 

1. ((Q — i£/)(n), r(n), (Q + iU)(n)) is an isotropic random field in <^-2©<^o© 
£,2- The fields <2(n) and U(n) are real-valued. 

2. (r(n),£'(n),5(n)) is an isotropic random field in <^o © <d0 © <^o with real- 
valued components. The components T(n) andB(n) are uncorrelated. The 
components E(n) and B(n) are uncorrelated. 

Proof. Let ((Q — if/)(n), T(n), (Q+iU)(n)) be an isotropic random field in <^_2© 
•sO © ^2, and let Q(n) and U (n) be real-valued. By Theorem 3 and reality condi- 
tions (23) and (30), we have E[ay^ w ] = for £ ^ 0, E[a±2j jn } = and 

E[ a T,em a T,£ r m'] = &U! '^mm'^T/, 
E[ a ±2,£m a ±Wm>) = ^U'^mm'^2,h ^4.) 
E[«r,fm«±2,£'m'] = 8g£f8 mm rCT,±2,£, 
E[a-2,£m«2,fm'] = ^1' $mm'C-2,2,t, 
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with 

oo CO 

£ (21 + l)C T/ + 2 £ (21 + l)C 2/ < oo. (35) 

Note that the second equation in (34) were proved for the general spin s by Geller 
and Marinucci (2008) in their Theorem 7.2. 

It is enough to prove that E[a E) £ m ] = E[a B) g m ] = and 

E[ax,g m aY,£'m'] = 8ee'8 m m'CxY,e (36) 

with 

£(2£+l)C^<°° (37) 

£ 

for all X, Y e {T,E,B}. Then, the second statement of the theorem follows from 
Theorem 3. 

The first condition trivially follows from (33). Condition (36) withX = Y = T 
is obvious. We prove condition (36) withX = Y = E. Indeed, by (33) and (30), 

^E/m^El'm'} = ^ (E[( a 2,£m + a -2,lm) i a 2,£'m' + a -2,£'m')} ) 

= ^[ fl 2/m^w] + E[a2/ ra a_ 2 /m'] 

+ E[a-2,£m«2,£'m'] + E [fl-2,lfn^w] ) 
= 2^ f £'^ mm '^ 2 ' e + Re ^-2,2,£)- 

Condition (36) with X = y = 5 can be proved similarly. 

Next, we prove condition (36) with X = T and y = 5. Indeed, 

Ef^r/m^B/m'] = ^(E[ar/m«2,fm'] — E[flr//n«-2,#m']) 

= 

by (30), which also proves that T(n) and B(n) are uncorrelated. Condition (36) 
for other cross-correlations can be proved similarly. 
Next, we prove (37) with X = E. Indeed, 

OO J oo 

£ (21 + l)C E>£ = - £ (21 + 1) (C v + ReC_ 2;V ) < °°. 

£=2 7 1=2 
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Condition (37) for X = B can be proved similarly. 

Next, we prove that E(p) is real- valued. It is enough to prove reality condition 

a>E,l -m = {-\) m aEi^. We have 

I, N 

-m — — 2 \ a 2,i ~m + a -2,£ -m) 

5(n) is real- valued by similar reasons. 

Finally, we prove that £'(n) and5(n) are uncorrelated. Indeed, E[.E (111)5(112)] = 
C£g(ni -112), because (r(n),£ , (n),5(n)) is an isotropic random field in <^o© <^o © 
£o- So, Ceb{(— ni) ■ (— n 2 )) = C£fi(ni ■ 112). On the other hand, 

C £5 ((-ni) ■ (-n 2 )) = E[£(-m)fi(-n 2 )] 
= E[ J E(n 1 )(-5(-n 2 ))] 
= -Cbb(iii -n 2 ), 

because E(—n{) =£ , (ni) and5(— ni) =— B(n{). Therefore, Ceb(^i • ^2) = 0. 

Conversely, let (7 1 (n) , E (n) , 5(n) ) be an isotropic random field in <^o © <^o © <^o 
with real- valued components, let the components T(n) and B(n) be uncorrelated, 
and let the components E(n) and B(n) be also uncorrelated. Solving system of 
equations (33), we obtain 

It is obvious that E[a± 2 #»] = 0. We have to prove (34), (35), and (30). The first 
equation in (34) is obvious. The second equation is proved as follows. 

E[«2/m«2/m'] = E[(— fl£ ! # n + a#/ m i) {—aE/'m 1 — «B/'m'i)] 
= <5«' <5 mm ' (C^ + Cg^) , 

because £(11) and 5(n) are uncorrelated. Proof for negative coefficients is similar. 
The third equation in (34) is proved as follows. 

E [flr>^/w] = E[ar ; £ m (— a E /' m i — «s/' m 'i)] 
= -8 u >d mm iC TE ,£, 
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because T(n) and B(n) are uncorrelated. Proof for negative coefficient is similar. 
The fourth equation in (34) is proved as follows. 

E[a-2/m«2,fm'] = E[( — fl£^ m _ ^B,lnS) ( _ «E,£W — a B,errX)] 

because E(n) and B(n) are uncorrelated. 
Because C21 = Cg £ + Cb e, we have 

00 00 00 00 

£ (2* + 1 )C T)i + 2 £ {li + l)C 2 ,e = £ (2£ + \)C T ,i + 2 £ {li + 1) (C £/ + C B/ ) 

e=o e=2 e=o t=i 

< 00 

which proves (35). The reality condition (30) is proved as 



a-2,£m = ~ a E,tm~^B,tm l 
= —aE,£m + UB,tnA 

= -(-l) m a Ej £ - m + (-\) m a B ,£ - m i 

— (—I) 1 "®!/ -m- 



□ 



In the so called Gaussian cosmological theories, the random field (r(n),£'(n),5(n)) 
is supposed to be Gaussian and isotropic with real-valued components. Let r/^oy, 
£ > 0, 1 < j < 3 and r)i m j, £> 1, 1 <m<£, 1 < j < 6 be independent standard 
normal random variables. Put 



C-&»j j 



7]£ 0j , m = 0, 

^(^m2j-l+^m2;i), m > 0, 



where £ > 0, < m < £, and 1 < j < 3. Now put 



1 /2 

O-TM = (Qr,i) C 



«£,&n — 1— 7TF7 bfrnl + ^El ~ Wm2j 
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for m > and ax£- m = (— l) m «xim for m < and X G {r,.E,.B}. The random 
fields 

f=0m=-l 
e=2m=-t 

satisfy all conditions of the second statement of Theorem 4. The random fields 
(Q±i£/)(n) can be reconstructed by (38), (28), and (29). By Theorem 4, ((£>- 
iU) (n) , T (n) , (<2 + iU) (n) ) is an isotropic Gaussian random field in <^_2 © £jo © £2 • 
The fields Q(n) and U (n) are real-valued. 

Finally, we note that Kamionkowski et al (1997) proposed a different for- 
malism for computations of the polarisation field on the whole sky. Instead of 
spin-weighted harmonics s Y^ m , they use tensor harmonics Yf and Y? which are 
related to the spin- weighted harmonics as follows. 

Yf m = (-2^ m e- © e_ + 2 Y( m e + © e+) , 

Y f m = t4f (-2*&»e_ © e_ - 2 Ye m e+ © e+) . 
iV2 

This formalism is also used by Cabella and Kamionkowski (2005), Challinor 
(2004, 2009), Challinor and Peiris (2009) among others. An excellent survey 
of different types of spherical harmonics may be found in Thorne (1980). 
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